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A subject  of  considerable  interest  to  the  C 


3 


community  is  the 


possible  effect  on  long-range  ELF  communication  of  widespread  ■ 

ionospheric  disturbances  of  nuclear  origin.  Predictive  ] \ 

• ’■  I 

calculations  at  the  present  time  are  based  on  the  so-called  | j 

"great  circle  WKB  approximation,"  in  which  the  properties  j j 


• • • j a 

of  the  entire  earth-ionosphere  waveguide  are  taken  to  be  those  , i 

that  exist  along  the  great  circle  between  transmitter  and  j ^ 

receiver,  and  the  properties  vary  slowly  in  the  horizontal  j j 

direction,  so  that  the  WKB  approximation  is  applicable.  j j 

4 \ 

\ 

1 

The  shortcomings  of  the  assumption  that  the  waveguide  is  1 

^ : 

laterally  homogeneous  have  been  recognized  for  a long  time.  ! ■ 

Because  of  the  mathematical  difficulties  involved  in  allow-  i ; 

1 

ing  for  lateral  variations,  the...;  have  been  very  few  theoret-  * 

ical  treatments  of  the  problem.  The  attempts  which  have  been  j l 

made  have  dealt  with  ionospheric  disturbances  of  limited  ) | 

i 1 

spatial  extent  [1,2]  which  are  amenable  to  scattering  and  dif-  l -i 

fraction  theory.  j j 


Here  we  are  concerned  instead  with  situations  involving  j 

transition  regions  where  conditions  vary  so  rapidly  that  a 
WKB  approximation  of  the  propagation  may  not  be  valid.  Such 
situations  occur  naturally  at  the  day-night  terminator  and  in 
the  polar  regions  during  large  solar  proton  events . In  this 
report,  we  investigate  the  effects  of  such  transition  regions  < 

I 

or  "edges"  on  ELF  propagation  along  and  through  them.  j 


T 


SECTION  2.  PHYSICAL  MODEL  FOR  NONUNIFORM  WAVEGUIDE 


It  is  well  known  that  waves  in  the  ELF  band  propagate  in  the 
earth- ionosphere  waveguide,  penetrating  to  ionospheric 
D-region  altitudes  and  higher,  depending  on  ionospheric 
conditions.  At  ELF,  only  the  lowest  (TEM)  waveguide  mode 
is  nonevanescent . The  mode  is  characterized  by  three 
parameters,  which  are  determined  primarily  by  the  ionization 
height  profile  of  the  ionosphere.  These  parameters  are  a 
relative  phase  velocity  c/v,  where  v is  the  phase  velocity  and 
c is  the  speed  of  light,  a loss  rate  a in  dB  per  thousand  kilo- 
meters of  propagation,  and  an  excitation  factor  AQ.  There  are 
a number  of  full-wave  numerical  methods  for  calculating  these 
parameters  when  the  ionosphere  is  assumed  to  be  vertically 
inhomogeneous  but  laterally  uniform.  There  is  no  full-wave 
numerical  method  for  determining  equivalent  parameters  in 
the  presence  of  lateral  inhomogeneities.  It  is  generally 
assumed  that,  in  a laterally  nonuniform  waveguide,  the 
propagation  may  be  characterized  by  local  parameters  with 
values  appropriate  to  a laterally  uniform  ionosphere  with  the 
local  vertical  profile.  This  assumption  forms  the  basis  for 
a number  of  approximate  solutions  for  propagation  in  a non- 
uniform  waveguide  [3].  Most  of  these  approximate  treatments 
require  that  the  lateral  variations  are  slow  enough  that  a 
WKB  approximation  may  be  used  to  calculate  the  lateral 
variation  of  the  ELF  fields.  Here,  we  will  be  interested  in 
propagation  in  the  vicinity  of  regions  where  this  criterion 
is  not  satisfied. 

The  problem  of  calculating  ELF  propagation  in  an  earth- 
ionosphere  waveguide  in  which  the  ionization  height  profile 
varies  along  the  propagation  path  is  very  difficult  mathe- 
matically. Even  in  the  case  of  horizontal  uniformity,  the 
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diffuseness  of  the  ionospheric  "boundary"  necessitates  a full 
wave  numerical  treatment  of  the  problem.  However,  once  such 
a numerical  solution  has  been  obtained,  it  is  possible  to 
replace  the  stratified  waveguide  by  a mathematically  equivalent 
sharply  bounded  waveguide  which  yields  the  same  values  of  the 
quantities  c/v,  a,  and  Aq.  This  is  effected  by  a proper  choice 
of  the  height  of  the  waveguide  and  the  surface  impedance  of 
its  upper  wall. 

In  principle,  it  should  also  be  possible  to  replace  a laterally 
varying  earth- ionosphere  waveguide  with  a mathematically 
equivalent  sharply  bounded  waveguide.  The  height  and  the 
surface  impedance  of  the  upper  wall  would  obviously  both  be 
functions  of  position.  In  order  to  keep  the  problem  mathe- 
matically simple,  we  will  consider  two  limiting  cases  within 
the  framework  of  such  a model.  In  the  first  case,  we  will 
assume  that  the  surface  impedance  of  the  waveguide  is  constant, 
but  that  the  height  of  the  guide  varies  in  the  lateral  direc- 
tion. In  the  second  case,  we  will  assume  that  the  height  is 
constant  and  that  the  surface  impedance  varies.  The  combina- 
tion of  the  two  cases  should  give  some  insight  into  the 
influence  of  lateral  gradients  on  ELF  propagation  in  the 
earth-ionosphere  waveguide. 


SECTION  3.  VARIABLE  HEIGHT  WAVEGUIDE  MODEL 


We  wish  to  model  a nuclear-disturbed  or  SPE-disturbed  earth- 
ionosphere  waveguide  in  which  the  waveguide  has  been  disturbed 
over  a very  large  area,  but  in  which  the  transition  between 
the  disturbed  and  undisturbed  portions  takes  place  in  a dis- 
tance small  compared  to  a wavelength.  As  discussed  in 
Section  2,  we  will  attempt  to  model  such  a situation  with  a 
sharply  bounded  waveguide  of  variable  height,  as  shown  in 
Figure  3-1.  For  mathematical  simplicity,  the  waveguide  is 
assumed  uniform  in  the  y-direction,  and  consists  of  three 
separate  sections.  The  outer  two  sections  are  semi-infinite 
uniform  waveguides  of  different  heights.  These  are  connected 
by  a central  section  with  a sloped  upper  wall.  The  waveguide 
of  greater  height  represents  the  undisturbed  portion  of  the 
earth-ionosphere  waveguide,  that  of  smaller  height  represents 
the  disturbed  portion,  and  the  central  section  represents  the 
transition  region.  The  lower  wall  is  assumed  to  have  infinite 
electrical  conductivity  and  the  upper  wall,  in  accordance  with 
the  discussion  of  the  previous  section,  to  have  constant  sur- 
face impedance. 

A cross  section  of  the  waveguide  is  shown  in  Figure  3-1, 
together  with  the  coordinate  system.  The  undisturbed,  tran- 
sition, and  disturbed  regions  are  designated  by  I,  II,  and  III, 
respectively.  The  origin  x = 0 of  the  transverse  coordinate 
is  taken  to  be  where  the  extrapolation  of  the  sloped  upper 
wall  intersects  the  lower  wall  in  Region  III.  The  source  is 
assumed  located  in  Region  I at  x = x , and  the  coordinates  of 

s> 

the  boundaries  of  the  transition  region  are  xQ  and  x^. 

f 

It  is  now  necessary  to  solve  the  wave  equation  in  each  of  the 
three  regions  and  then  to  apply  the  appropriate  boundary 
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conditions  at  x = x and  x = x,  . We  assume  the  source  is  an 
electric  dipole  antenna,  oriented  either  vertically  or 
Horizontally.  The  first  configuration  would  correspond  to 
a natural  source  of  ELF  radiation,  such  as  a lightning 
discharge,  and  the  other  to  an  articifial  source,  such  as 
the  Wisconsin  Test  Facility.  We  assume  that  the  source  is 
sufficiently  far  from  the  transition  region  that  only  the 
nonevanescent  TEM  mode  is  incident  on  the  transition  region. 
We  also  assume  that  no  other  modes  are  generated  in  the  tran- 
sition region  or  at  its  boundaries. 

In  the  two  uniform  regions,  the  TEM  fields  can  be  derived 
from  an  electric  Hertz  vector  which  has  only  a vertical 
component  11^.  The  electromagnetic  field  is  obtained  from  the 
Hertz  vector  through  the  equations 


£ = $ + k2  n i 

\3y  / y y 


(3-la) 


3 = 


i w 


A 


(3-lb) 


where  k = w/c  is  the  free  space  wave  number  and  Py  is  a unit 
vector  in  the  y-direction.  In  Region  I,  the  primary  excitations 
for  the  different  dipoles  are  given  by 


n 


VD 


= r—  cos  (kC  y)  H 


(1) 


kS„  ( (x- 


> 2 , 2 
x ) + z 


I J 

! 

i i 
i 1 
i ? 


j 


) i 


> 

: i 

1 i 


( 3-2a) 


where  VO,  HX  and  HZ  refer  to  a vertical  dipole,  a horizontal 
dipole  oriented  along  the  x-direction,  and  a horizontal  dipole 
oriented  along  the  z-direction,  respectively.  In  these 
expressions,  A0  is  a factor  depending  on  the  source  strength, 
hQ  is  the  height  of  the  guide  in  Region  I,  and  SQ  and  CQ  are 
the  complex  sine  and  cosine  for  the  TEM  mode  in  a uniform 
waveguide  with  the  properties  of  Region  I. 

Since  the  properties  of  the  waveguide  are  independent  of  z, 
it  is  useful  to  carry  out  a Fourier  transformation  in  this 
variable.  Thus,  we  write 


ny (x,y,Sz) 


-ikS_  z 


z)  e 


dz 


(3-3) 


with  similar  expressions  for  the  components  of  E and  B.  The 
transforms  of  the  primary  excitations  given  by  Equation  (3-2) 
can  be  written  in  the  form 


ny  (x,y,Sz) 


- A JDIP 


(S  ,S  ) cos  (kC  y)  e 

U c v 


2 2 

-ik(S^-S^)  (x-xs) 


(3-4) 


DIP 

where  the  function  F corresponding  to  the  different  dipoles 
is 


FVD  = (8*  - S*> 
o z 


-1/2 


( 3-5a) 


FHX  = 1 


,HZ 


s (s2n 

z o 


-1/2 


(3-5b) 

(3-5c) 
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In  Region  I,  the  total  field  consists  of  the  primary  excitation 
radiated  by  the  source  plus  the  wave  reflected  from  the  transition 
region.  Thus,  the  total  Hertz  vector  in  Region  I can  be 
written  as 


n.  = 


2 2 

Ao  DIP  ■ lk(So‘Sz)  (xs"Xl) 

jp  cos  (kCQy)  e ° 2 s x 


j 5 1/2  2 o 1/2 

-ik(S^-S2)  (x-x.)  ik(S^-S2)  (x-x.) 

e 02  1 + R e ° 2 1 


(3-6) 


where  R is  a reflection  coefficient  to  be  determined.  In 
Region  III,  on  the  other  hand,  there  is  only  an  outgoing  wave, 
so  that  the  Hertz  vector  can  be  written  as 


r 1 1 1 


A1  DIP  ik<So"Sj> 

~ FUiF  cos  (kCxy)  e ~ “ 


1/2 


(xs-xi}  -ik(S1-Sz) 
e 


1/2 


(x-x, 


o 


(3-7) 


where  and  C-^  are  the  complex  sine  and  cosine  appropriate  to 
Region  III  and  A^  is  a constant  to  be  determined. 

A relationship  exists  between  and  Sq  as  a consequence  of  the 
assumption  that  the  upper  wall  has  a uniform  surface  impedance. 
From  Equation  (3-la)  and  (3-lb) , the  tangential  field  components 


«n **■  fwDtifi.  annptinnNMmp*  » ukrtw^ish!*.,  -» * .o^****  * mtoinum  turn  * mmmMI  *i«fi  #4§&Plfc|MP***9»*l^R 
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It  now  remains  to  construct  the  solutions  in  the  transition 
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It  can  be  expected,  for  the  kinds  of  disturbances  with  which 
we  are  concerned,  that  the  change  in  height  across  the 
transition  region  is  perhaps  a few  tens  of  kilometers  while 
the  lateral  distance  over  which  it  takes  place  is  of  the  order 
hundreds  of  kilometers.  Thus,  from  Figure  3-1, 


— <<  1 


(3-17) 


and  the  angle  <j>  is  restricted  to  the  range  0 £ 4>  £ <j>  <<  1. 

Furthermore,  we  know  that  at  the  ground  (<$>  = 0)  the  tangential 
components  E^  and  of  the  electric  field  must  vanish.  Thus, 
the  power  series  expansions  for  u and  w must  have  the  forms 


u = G<£)  - g(K)  . 


w = F ( 0 4>  - f (?)  f-  + . . • 


( 3-18a) 


(3-18b) 


When  these  solutions  are  substituted  into  Equation  (3-14) , we 
obtain  the  differential  equations 


"1  1 2 1 
G + | G + (1  - S^)  G - ~ g - 0 


" 1 ' 2 1 
F +iF  +(i_s^)F_i_f=0 


(3-19a) 


(3-19b) 


In  order  to  obtain  a simple  differential  equation  for  one  of 
the  four  functions  appearing  in  these  two  equations,  it  is 
necessary  to  obtain  two  additional  relations  among  them.  These 
are  provided  from  the  known  surface  impedance  of  the  upper  wall. 
The  components  of  the  electric  and  magnetic  fields  in  Region  II, 
as  obtained  from  Equations  (3-15),  (3-16),  and  (3-18),  are 


*r  - - k 


u 3 + s* F /♦ 


(3-20a) 


Ef  = - k2 

<P 


(G'  + slxF') 


EZI  = 1 ^ SZ  (1  “ S2)  F * 


B 


II 


B 


P 

II 


■ £ 8*(°'  + k p) 


2 

- i ~ a - s“)  g 

GJ  z 


The  surface  impedance  at  the  upper  wall  (<J>  = 


V is 


1 

o 


E 


ii 


,n 


CB 


II 


cB 


II 


This  leads  to  the  two  relations 


(1-S2)F 


(G  +FA  ) 


. A..  = c2 e 

<f>  n -o'l 


(g/5+S‘F  ) 


(1— Sg ) G 


_1_  = c2£ 

<t>  n co^i 


(3-20b) 


(3-20c) 


(3-21a) 

(3-21b) 


(3-22) 


(3-23a) 


(3-23b) 


where  the  last  equality  in  each  equation  follows  from  Equation 
(3-10)  and  the  geometry  of  Figure  3-1. 
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In  order  to  apply  the  boundary  conditions,  we  must  evaluate 
Equations  (3-26)  and  (3-28)  at  £ = The  results  are 


G (C1) 


G-  (£ , *A 

^ ' 5"  2 2 

Vo  <S1-Sz> 


(3-29a) 


G(CX) 


G(^0) 


1 + 


G Eq) 

GUP“ 


° <srsz> 


h 

in  r~ 
hl 


(3-29b) 


The  solution  is  completed  by  making  the  fields  continuous  at 
x = xq  and  x = x,  . Since  the  slope  of  the  upper  wall  in 
Region  II  is  very  small,  the  p-component  of  the  fields  is  nearly 
horizontal.  Thus,  we  require  that 


B 


B 


II 

P 

II 

z 


BxII(xo> 

BzII(x0> 


(3-30a) 


II 

P 

II 

z 


(Pj_) 

(Px) 


Exlxl> 

BX> 


( 3-30b) 


The  first  of  these  equations  leads  to  the  relation 


G (50) 

GUq} 


i (s: 


1/2 


(3-31) 


where  the  left-hand  side  of  Equation  (3-30a)  has  been  calculated 
from  Equations  (3-21a)  and  (3-21b),  making  use  of  Equation  (3-23), 
and  the  right-hand  side  has  been  calculated  from  Equation  (3-8) 
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er 
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* 


with  II*11  given  by  Equation  (3-7)  . Carrying  out  the  same 
procedure  for  Equation  (3- 2 3b) , we  obtain 


& 

& 


G (?1) 

gT^T 


= - i s2  - s2)1/2  a + R< 
1 So  *V  (I-  R) 


(3-32) 


The  last  equation  allows  us  to  calculate  the  reflection 
coefficient 


R = 


G (C1)/G(51)  + i(S~-S^) 
g'(?1)/G(C1)  - i(S^-S \) 


1/2 


(3-33) 


which,  making  use  of  Equations  (3-28),  (3-29),  and  (3-31), 
becomes 

1/2 


9 2 9 2 2 2 

sfh,  (s*-sn  - s^i  (sr-sj) 
R - 1 1 o Z O O.  - 1 z 


1/2 


(i-iy) 


2 2 2 


2 .2 


s^h,  (s>s,)  + s"h„(st-s:) 

1 i o z 


1/2 


o o'  1 ZJ 


(1-iy ) 


(3-34a) 


r* 


where 


Y = 


EoSl 


2 2 

(srv 


1/2 


Jin 


S Jin 


vi 

hisi 


( 3-34b) 


Finally,  we  need  a relation  between  A-^  and  Aq,  the  amplitudes 
of  Hertz  vectors  in  the  disturbed  and  undisturbed  regions, 
respectively.  This  is  obtained  by  making  B continuous  at 
x = xQ  and  x = x^.  These  two  conditions  may  be  expressed  as 


Bz  <*1> 


(3-35) 
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Carrying  out  this  operation  and  making  use  of  the  various 
relations  above,  we  obtain 
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2 2 
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(3-36) 


All  quantities  necessary  for  the  calculation  of  the  electromag- 
netic field  in  Region  I and  II  have  now  been  obtained.  In 
practice,  the  fields  are  detected  either  by  whip  antennas, 
which  measure  the  vertical  electric  field,  or  by  loop  antennas, 
which  measure  the  horizontal  magnetic  field.  To  obtain  numerical 
values  for  these  quantities,  it  is  necessary  to  invert  the 
appropriate  Fourier  transforms.  These  results  will  be  presented 
in  a subsequent  section. 
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SECTION  4.  VARIABLE  IMPEDANCE  WAVEGUIDE  MODEL 
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As  in  Section  3,  we  again  assume  a waveguide  which  consists 
of  three  separate  sections  which  are  uniform  in  the  y-direction. 
However,  we  now  assume  that  the  height  of  the  waveguide  is 
constant,  but  that  the  surface  impedance  is  not.  The  outer 
two  sections  are  again  semi-infinite  uniform  waveguides, 
with  different  surface  impedances.  They  are  connected  by  a 
central  section  of  length  L whose  surface  impedance  varies 
continuously  between  the  two  outer  values. 


The  Hertz  vectors  in  Region  I and  Region  III  are  again  given 
by  Equation  (3-6)  and  (3-7) , respectively.  However,  there  is 
no  longer  a functional  relationship,  such  as  Equation  (3-12) , 
between  the  eigenvalues  in  the  two  regions.  The  quantities 
S and  S are  independent  and  considered  specified.  The  Hertz 
vector  in  the  transition  region  satisfies  the  wave  equation 


with 


sV1 


3x 


+ k 


S2(x) 


(xQ  < X < Xf) 


(4-1) 


S(xo) 


1 


S (x1) 


S 


o 


(4-la) 


In  analogy  with  the  sloped  waveguide  of  the  preceding  section, 
we  assume  the  impedance  of  the  transition  region  varies 
linearly  in  the  transverse  direction,  i.e.. 
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(4-2) 


7 7 7 7 

s2(x)  = s2  + <s2  - s2) 


With  S (x)  given  by  Equation  (4-2),  the  solution  of  Equation 
(4-1)  is 


u 


.II 


- B u1/3  [j1/3(u)  + 3 J-i/3<u>] 


where 


u = (ate) 


k = k [s2(x)  - S2  ] 


1 1/2 


a « ? 


3 k2(S2-S3) 


and  B and  3 are  two  constants  to  be  determined. 


(4-3) 


(4-4a) 

(4-4b) 

(4-4c) 
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The  boundary  conditions  are  again  that  the  tangential  electric 
and  magnetic  fields  must  be  continuous  at  the  boundaries  of 
the  transition  region.  This  requires  IIy  and  to  be  contin- 
uous at  x = xQ  and  x = x^.  From  the  continuity  of  the  ratio 
of  these  quantities  at  x * xQ,  we  obtain 


_ i J1/3<V  ♦ J-2/3(uo) 
J2/3  <uo>  ' i J-l/3(uo’ 


(4-5) 


and  from  the  continuity  of  the  same  quantity  at  x = 

[Jl/3(ul)~1  J-2/3(ul)]  + 3 [ J-l/3  (ul)-f:L  J2/3(ul)] 
[J1/3(U1)+1  J-2/3(ul)]  + B j^-l/S^l5"1  J2/3(U1}]  (4-6) 
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Finally,  from  the  continuity  of  n at  x = xQ  and  x - x^, 


A!  ° B u 


AQ (1+R)  = B 


‘o/3  [Jl/3<“o>  + 6 J-l/3(uo>; 

(4-7a) 
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from  which 
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(4-8) 


With  the  above  expressions,  the  Fourier  transforms  of  the  fields 
may  now  be  inverted  in  a straightforward  manner.  However,  in 
view  of  the  idealization  of  the  model,  it  seems  more  reasonable 
to  look  at  a limiting  case  where  considerable  simplification 
occurs  in  the  complexity  of  these  expressions  and,  therefore, 
also  in  the  necessary  computational  labor.  This  is  the  limit 
in  which  the  lateral  extent  of  the  transition  region  becomes 
small,  i.e.,  L ■+  0.  In  this  case,  the  Bessel  functions  in  the 
above  expressions  may  be  approximated  by  their  small  argument 
expansions,  and  Equations  (4-6)  and  (4-8)  become 
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Numerical  results  based  on  these  simpler  expressions  are 
presented  in  the  next  section. 


SECTION  5.  NUMERICAL  RESULTS  AND  CONCLUSIONS 


The  results  of  the  preceding  two  sections  have  been  used  to 
obtain  numerical  values  for  the  magnetic  field  along  various 
great  circle  paths  through  the  transmitter,  for  the  ionospheric 
parameters  listed  in  Table  1.  These  results  are  compared  in 
Figures  5-1  to  5-6  with  those  calculated  from  the  great  circle 
WKB  approximation  for  the  same  paths.  The  great  circle  WKB 
results  are  obtained  from  the  following  expressions: 

Region  I 


VD  = 1/2  FI 

bwkb  So 


(5-la) 


S0/2  F1  cos  0 


(5-lb) 


S2/2  F1  sin  0 


(5-lc) 


where 
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and  r is  the  distance  measured  along  the  propagation  path. 
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Region  III 


VD  _ 3/4  -1/4  III 

BWKB  " So  S1  F 


( 5-3a) 


„HX  _ / 4 _-l/ 4 pIII  n 

bwkb  - So  si  F cos  0 


( 5-3b) 


aHZ  = s7/4  s“1/4  F111  sin  0 

bwkb  so  si  F sin  0 


( 5-3c) 
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Table  5-1.  Ionospheric  Parameters  for  the  Model  Waveguides 
Used  in  the  Calculations* 
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Figure  5-6.  Mag 
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In  these  expressions,  6 is  the  angle  of  the  great  circle  path 
with  respect  to  the  x-axis.  Thus,  0 = 0 is  the  path  perpen- 
dicular to  the  transition  region  and  0 = 90°  is  the  path 
(entirely  in  Region  I)  parallel  to  the  transition  region. 


In  the  figures,  we  have  plotted  the  quantity 

.DIP 


20  log10 


B 


B 


DIP 

’WKB 


as  a function  of  distance  from  the  transmitter,  for  a trans- 
mitter frequency  of  50  Hz.  Figures  5-1  to  5-3  show  the  results 
for  the  variable  height  waveguide,  the  parameters  for  which  are 
given  in  Table  1.  Figures  5-4  to  5-6  present  the  results  for 
the  variable  impedance  waveguide , the  parameters  for  which  are 
also  given  in  Table  1.  The  models  have  been  constructed  so 
that  the  phase  speeds  and  attenuation  rates  are  the  same  for 
each  in  Regions  I and  III.  The  modal  parameters  and  height  of 
the  undisturbed  region  are  representative  of  an  ambient  day- 
time ionosphere.  The  parameters  for  the  disturbed  region  are 
representative  of  a nuclear  disturbance  with  an  ionization 
intensity  factor  f = 10  . This  factor  is  defined  as 


, (FY) 

' rV'2 


(5-6) 


30 


where  FY  is  the  total  fission  yield  in  megatons,  t is  the 
time  in  seconds  after  the  burst,  and  R is  the  radius  in 
kilometers  over  which  the  debris  is  assumed  to  be  uniformly 
spread.  In  this  report,  our  primary  concern  is  not  with 
the  modelling  of  a particular  disturbance,  but  with  the 
accuracy  of  various  methods  of  carrying  out  propagation  cal- 
culations in  the  presence  of  a sharp  lateral  inhomogeneity. 

It  can  be  seen  from  the  figures  that  the  results  for  the  two 
models  are  not  very  different  from  each  other,  and  that  both 
are  In  agreement  with  the  WKB  results  to  within  one  or  two 
dB  at  distances  of  several  thousand  kilometers.  Since  the 
total  attenuation  in  this  distance  is  several  dB,  we 
conclude  that  the  great  circle  WKB  approximation  is  quite 
adequate  in  the  calculation  of  ELF  propagation  through  wide- 
spread disturbances,  even  in  the  presence  of  "sharp" 
boundaries . 

Of  cou..  , it  must  be  realized  that  the  ionospheric  models 
employee  n the  calculations  have  been  rather  artificial, 
since  a real  ionospheric  disturbance  would  change  both  the 
height  and  '.mpedance  of  the  upper  boundary.  It  is  possible 
that  calculations  based  on  a more  realistic  model  would 
reveal  larger  discrepancies.  This  remains  a matter  for 
further  investigation. 
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